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ABSTRACT 

Temperature  tables  cure  given  for  the  calculation  of  tempera¬ 
tures  in  a  two-layer  slab  as  a  function  of  distance  from  one  face 
and  time.  Heat  is  transferred  at  the  interface  and  at  the  two 
outside  faces.  The  complete  solution  is  a  linear  combination,  V, 
and  Vg,  tabulated  in  Volume  1  and  Volume  2,  respectively. 
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INTRODUCTIOR 

The  problem  considered  here  is  determining  the  temperature 
in  a  two-layer  slab  with  heat  transfer  at  both  outside  faces 
and  a  "contact  resistance"  at  the  interface.  The  initial  tem¬ 
peratures  of  the  two  layers  are  constant  and  the  outside  stag¬ 
nation  temperatures  are  constants.  These  are  special  cases  of 
the  problems  of  references  1  and  2.  Details  of  the  solutions 
are  given  there. 

THE  PROBLEM 

The  equations  to  be  satisfied  are 

U1XX  "  klUlt  Tor  0  <  X  <  V  t  >  °>  (1) 

U2XX  *  k2U2t  for  h  <  X<  V  *  >  °>  (®> 

U1  -  ua  for  C  <  X  <  1^,  t  -  0,  (3) 

Ug  -  Ub  for  ^  <  X  <  Lg,  t  -  0,  (*) 

K1U1X  *  'hi  *Ui  "  V  for  X  *  0*  t  >  0,  (5) 

K1U1X  "  K2U2X  “^c  (U1  "  V  f°r  x  *  L!*  t  >  °>  (6) 

K2U2X  *  -ho  (U2  '  V  for  x  “  V  t  >  0.  (7) 

The  temperatures  in  degrees  are  given  by  Eq.  (17).  V^X^t) 
is  given  in  Table  1.  Vg(X,t)  is  given  in  Table  2.  Vj(X,t)  is 
the  value  of  Vg  in  Table  2  for  the  parameter  and  variable  cor¬ 
respondence  of  Eq.  18. 

is  tabulated  in  this  volume,  and  Vg  is  tabulated  in  Volume  2. 
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is  directly  the  temperature  for 

ui  =  1,  ua  =  0,  =  0,  Uo  =  0  . 

is  directly  the  temperature  for 

Ui  -  °*  ua  =  1,  ub  =  0,  Uo  =  0  . 


is  directly  the  temperature  for 

U  =  0,  U  =  0,  U  =  1,  U  -  0  . 
i  a  b  ’  o 


THE  SOLUTION 

Changing  variables  to  x,r  Eqs.  (l)  through  (7)  can  be  written 


'lxx  -  V  for0<x<o1. 

,  r>  o, 

(8) 

.  2,.  for  a.  <  x  < 

2xx  "  **2  <>r  *■ 

l,r>  0, 

(9) 

Ul  "  for  0  <  x  <  0^, 

r«  o, 

(10) 

Ug  «  for  <  x  <  1, 

r«  o. 

(11) 

'lx  "  'ai(Ui  -  V  for  x  ■ 

o,  r>  o, 

(12) 
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Ulx  "  °3U2x  "  'ac^Ul  •  V  forx-Vr>  °’ 

U, 


2x 


-ao(U2  -  UQ)  for  x  ■  1,  r>  o. 


(13) 

(1M 


It  can  be  seen  that  Eqs.  (8)  through  (14)  will  be  satisfied  by 

u(x,t)  =  UjV^x.r)  +  uav2(x,  r)  +  utv5(x,r)  +  uvjx,r)  (15) 

'u1(x,r)  0  <  X  < 

Ug(x,r)  <  x  <  1 

and  V, ,  V_,  V,,  V,.  are  the  solutions  of  Cases  I  -  IV: 

1  2  3  **■ 


where  U(x,f)  = 


(16) 


Case  I 
Case  II 
Case  III 
Case  IV 


U  -  1,  U  •  0,  U  =  0,  U  *  0, 
i  ’a  ’  b  o 

U  -  0,  U  ■  1,  U  «  0,  U  «  0, 

1  .  8  D  O 


u  -  o,  u  *  o,  u  «  1,  u  -  0, 

i  a  b  o 


ui  »  °,  ua  -  «  0,  Uq  •  1. 


It  can  also  be  seen  that 

u(x,t)  -  uo  -  (ut  -  uo)  vx(x,t)  +  (ua  -  uo)  v2(x,r) 
+  (ub  -  uo)  v3(x,r), 

and  that 


(17) 


(18) 


This  symmetry  is,  with  a  few  exceptions,  preserved  in  construc¬ 
ting  tables  of  W1  and  Vg;  i.e.,  for  many  values  of  aQ,  0^,  a,,  0^, 

Og,  (Xj  found  in  the  tables  the  "symmetric- Image”,  0^,  aQ,  aJay 

. ..,  may  also  be  found  in  the  tables.  Values  of  V^  are  to  be  found 

from  corresponding  entries  in  Vg  as  given  by  Eq.  (18).  In  some 

cases  the  desired  entries  are  outside  the  range  given  but  many  of 
them  are  in  regions  where  they  will  be  closely  given  by  the  last 
table  entry. 
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Taking  the  Laplace  transform  of  Eqs.  (8),  (9),  (12),  ( 13)  and 
(l4),  we  have 


CD 


U 


lLxx 


U2Lxx  =  ~a2Y>2v2L  "  °2 


(s  =  -b2,  U1L  s  J  e-sr"  U^XjfJdT) 

=  -b2U1L  -  U&  for  0  <  x  <  a^, 

for  <  x  <  1, 

b2uiLx  =  ai(Ui  +  b2uiL)  f°r  X  “  °» 

UlLx  =  a3U2Lx  =  'ac(UlL  -  U2L>  f°r  X  “  °L» 


U2Lx  “  ao(Uo  "  U2L)  for  X  "  1- 


The  solution  to  ( 19)  and  (20)  is  of  the  form 
U1L  =  c1  sin  bx  +  c2  cos  bx  +  P1(x,b) 

and 


(19) 

(20) 
(21) 
(22) 

(23) 

(24) 


U2L  *  c^  sin  Ogbx  +  c^  cos  Ogbx  +  P2(x,b).  (25) 

Substitution  in  Eqs.  (21)  through  (23)  gives 

-  o^Cg  -  aiui  +  aib2P1(0,b)  -  bS^to.b)  -  A1(b)  (26) 

c^b  cos  c^b  -  Cgb  sin  C^b  -  c^cyXjb  cos  ayXjb  +  c^agOtjb  sin  OgC^b 

-  *Plx(ai'b)  +  a3P2x(aL,b)  £  A2(b)  (27) 


c1(b  cos  a^b  +  ac  sin  Ojb)  +  Cg(-b  sin  a,  cos  atjb)  -  c^ac  sin  OjOgb 
-c^ac  cos  a^atgb  «  0^2(0^^)  -  0^(0^  ,b)  -  plx(°t1»b)  *  Aj(b) 

(28) 
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Cj( CKgb  cos  a^h  +  Oo  sin  a2b)  +  c^-a^b  sin  atgb  +  ctQ  cos  Ofgb)  * 

(29) 


a  u 

o  o 


-0oP2(i,b)  -  p^d.b)  -  -V  *  Vb)- 

D 

The  solutions  to  Eq.  (26)  through  (29)  will  be  of  the  form 
N  .(b) 

ci  =  DTbT  for  1  “  2»  3>  4- 

The  mcst  convenient  form  for  writing  N  (b)  and  D(b)  is  found  by 
using  the  two-row  two-column  Laplace  minors. 

Define 


'U 


ai,J 

al,j+l 

’t>ro 

ii 

ai,j  ai,j+l 

ai+l,J 

ai+l,J+l 

ij 

ai+2,j  ai+2, j+1 

D(b)  = 


all 

**12 

0 

0 

a21 

a22 

a25 

*2k 

a51 

a52 

a55 

a34 

0 

0 

a43 

a44 

Then 


D(b)  =  A^"  -  A^  A^ 

v  '  11  *33  11  23 


(50) 


N^(b)  sin  bx  +  N^(b)  cox  bx  = 

A55  [A1  a22  sin  1t>x  *  a21  cos  bx] 

+  A2 

|a11  cos  bx  -  a12  sin  bxj 

}  4[A25  Ai 

a^  cos  bx  -  sin  bxj 

A5 

a12  sin  bx  -  a^  cos  bxj 

}  4«a‘ 

cos  bx  -  a12  sin  bxj 

(51) 
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N^b)  sin  a2bx  +  N1+(b)  cos  d^bx  *  A2i{Ai[aU4  sin  a2bx 


-%3  cos  a2bx 


]}  *  4f*2[- 


a^3  cos  a2bx  -  &kh  sin 


<¥“] 


v 

L  a1  L 

- 

a24  sin  Qtgbx 

-  a25  cos  a2bx| 

j  +  *ii(A 

3  aU4 

■a45  cos  (Sgbxj  +  Ajjjajj  cos  Ogbx  -  sin  Ogbxjj 


The  solutions  are 


^{ajuajl-c^)]  ♦a0a,(i«y»e)  -<W»C*) 
1  ’  '  ai«c  +  ao[aiacU-al4<xla3)  +  a3(ai  +  V] 

a)  lb  e  _bn 

Jl  [Vbn>  8l”  V  *  VV  008  Vj-TTiT 


-  2 


U0(x,f) 


_ VolK  -  Vl 

ai°c  *  ao pi0^  1_al4ala3^  * 


-2  »Il[" 


5(bQ)  sin  c^bnx  +  N^b^  cos  a^x 


-b2r 
be  n 


(32) 


(33) 


(34) 


where  the  bR' s  are  the  positive  ordered  zeroes  of  D(b)  and 

D'(b  )  dD(b)  at  b  ■  b  .  Further  details  are  given  in  Ref.  1  or 
db 

Ref.  2.  The  functions  actually  tabulated  are  the  and  V2of 

Case  I  and  Case  II,  respectively.  Equation (l8) can  be  used  to  find 
and  actual  t< 

solutions  were: 

For  V^,  ■  0,  Pg  ■  0  and  thus 

A1  *  ai’  A2  "  °’  A3  *  °»  kh  “  0 


Vj  and  actual  temperatures  are  found  by  using  Eq.  (17).  The  particu- 
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For  Vg,  P1(x,b)  *  -^2  (1  -  cos  bx),  P2  =  0 
b 


and  thus 


sin 


.  .  .  **1  .  8inbai  -  1-cos  bx 

Ax  «  0,  A2(b)  «  - 5 -  ,  A  (b)  =  - £ -  +  ac  - g - » 


\  -  o 


Characteristic  numbers  (bn,  n«l,2,  ...  ,  the  positive  zeroes  of  D(b)) 

By  introducing  "variable  phase  angles"  into  the  function 

^♦(b)  *  ,  these  angles  being  monotone  increasing  functions  of 

b 

b,  we  can  see  that  D(b)  *  0  if  and  only  if  h^m.fb))  *  K[h2(m2(b))+b], 
K  *  pos.  const.,  where  the  m^  are  monotone  increasing  functions  of  b 
and,  the  h^  are — up  to  constant  multipliers — "oppositely  directed" 
tangent  functions.  This  yields  two  sets  of  b-values,  (A.)i-0,1,2, . . . 
and  (Bj)j^>,l,2, . . . ,  the  points  where  the  vertical  asymptotes  of  h^ 
and  h2  cut  the  b-axis,  with  the  property  that  between  two  adjacent 
values  there  lies  one  and  only  one  zero  of  D(b). 

The  only  caution  to  be  observed  is  to  note  that,  if  A^  ■  B^ 

(for  some  i,j)  then  the  zero  "between"  A^  and  Bj  is  A^  «  Bj. 

The  location  of  the  (A  )  and  (B.)  reduces  to  the  location  of 

X  J 

the  S£gogg  of  functions  of  the  form  f(b)  ■  c^  cos(K'b)-c2b  sin(K,b), 
c^  non-neg.  consts.,  K'  a  pos.  const.,  and  the  zeroes  of  this 
function  occur  at  precisely  those  values  of  b  for  which 

ctg(k’b)  ■  K"b,  K"  a  non-neg.  const. 

Again,  we  have  two  "oppositely-directed"  functions,  only  one 
of  which  is  tangent  type,  so  we  can  find  Just  one  set  (Ci)i^5,l,2, . . . 

of  b-valuer  such  that  "between"  two  adjacent  ones  there  lies  one  and 
only  one  zero  of  f(b). 

A  similar  caution  must  be  taken  here — for  certain  values  of  a, 
the  zero  "between"  and  C1+1  can  actually  be  (or  C1+1).  In 

working  out  the  explicit  details,  the  cases  where  this  happens  are 
simple  to  observe. 


7 


NAVORD  REPORT  5362 
Part  7 


THE  TABLES 


The  values  of  the  parameters  and  variables  used  in  the  tables . 

A.  For  all  parameter-settings,  the  following  values  of  time, 
? ,  and  proportional  distance,  p,  are  used: 

r*  10  10  ,  10  ,  1,  10,  10  ,  HT,  10  , 

occurring  across  the  top  of  the  body  of  any  one  column  as 

-21  -11  1  11  21  31  51  51; 

•p  -0,  .2,  .4,  .6,  .8,  1.0;  0,  .2,  .4,  .6,  .8,  1.0. 

1st  layer  2nd  layer 

Distance,  x,  is  connected  to  p  by: 


x  *  pOL ,  in  1st  layer, 

1  (35) 

x  *  +  p( 1-0^) ,  in  2nd  layer 

We  might  note,  in  particular,  that  U(p^,f)  *  U(l^,T)  is  the 

answer  at  the  right  edge  of  layer  one  (p^  means  "a  p  from  the 
first  string  of  p's";  p'2'  means  "a  p  from  the  second  string  of  p's"), 

whereas  U(p^,  r)  ■  U (0^,f)  is  the  answer  at  the  left  edge  of 
layer  two;  x  itself  is  (X^,  in  both  cases. 

B.  Two  sets  of  a's  are  used: 


and 


0 


io:3. 

io"2, 

io'1. 

1, 

10 

10‘3, 

io'2, 

io'1, 

1, 

10 

io' 3, 

CM 

I 

O 

H 

IO'1, 

1, 

10 

(l  value) 
(5  values) 

io2,  io5  (7  value8) 

102,  103  (7  values) 


a,  10(-|)  (3  values) 

a2  a2  a2 

-  .01,  .05,  .1,  .2,  .5,  .8,  .9,  .95,  .99  (9  values) 
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2 

The  a2~f  the  cty  and  the  ct^-  set  stay  the  same.  One  column 
(for  one  case)  is  taken  up  by  a  sweep  of  the  nine  0^-  values.  Thus, 
the  sweep  of  this  set  occupies  the  first  735  columns  of  each  table. 


a0 

=  10'2, 

io"1. 

1, 

10 

ai 

O 

O 

II 

‘2,  10' 

-1 

* 

1,  10 

a 

c 

OJ  " 

1 

o 

H 

II 

10'\ 

1, 

10 

p 

The  sweep  of  these  (together  with  the  a y,  ay, 

above)  takes  1,680  columns  or  c.  736  through  c.2,4l5 
column  of  a  table). 


(4  values) 

(5  values) 

(4  values) 

o^-  sets 
(the  last 


The  order  of  sweep  is  the  inverse  of  the  ordering  above;  i.e., 

2 

the  ay  set  is  swept  out,  then  the  ay  set,  the  ac-set,  etc. 


The  accuracy — Symbols  A,  B,  C,  X,  U. 

(  q\ 

In  any  one  clump  of  answers  (a  sweep  of  pVP'  and  tr),  order 
the  answers  from  upper  left  to  lower  right,  in  the  natural  way, 
(i.e.,  sweep  the  first  row  from  left  to  right,  then  the  second  row, 
etc . ) . 


The  accuracy  of  the  first  answer  of  the  clump  is  denoted  by 
one  of  the  above  symbols.  This  symbol  also  indicates  the  accuracy 
of  following  answers  until  the  accuracy  changes,  when  a  new  symbol 
is  met. 

The  (symbol-accuracy)  connection  is  meant  to  be: 

A  -  All  4  digits  have  meaning  (4-digit  accuracy) 

B  -  Right-most  digit  is  meaningless  (3-digit  accuracy) 

C  -  The  2  right-most  digits  are  meaningless 
X  -  The  3  right-most  digits  are  meaningless 
U  -  Your  guess  -If  you  guess  between  0  and  1  - 
as  good  as  the  one  printed! 

(The  right-most  meaningful  digit  may  be  wrong,  but  by, 
at  most,  1  unit). 

Unfortunately,  there  is  at  least  one  section  of  the  tables 
in  which  the  above  table  is  not  to  be  trusted.  It  is  c.2l40. 

Case  II,  with  oq  ■  10  •,  -  0.01,  a,  *  1,  a|  « 1000.  ,  at ^  »  .0001, 

0^  -  .99. 
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The  loss  of  meaning  of  the  symbol  (here.  A)  or,  equivalently, 
the  loss  of  accuracy  of  the  table-entries,  was  due  to  the  oversight 
of  an  obvious  source  of  round-off  error  on  the  part  of  the  programmer 
and  was  Initially  confined  to  Isolated  clunps  In  the  tables.  However, 
many- -If  not  all— of  the  worst  clumps  were  re-done,  after  the  program 
had  been  corrected. 

These  errors,  when  they  are  obvious,  are  made  obvious  by  a 
certain  local  monotonicity  (In  either  p  or  t)  being  violated  (in 
the  case  of  c.  21 40,  for  fixed  and  smallest  value  of  r,  in  the 
second  layer,  the  table-entries  vhould  be  zero,  clearly,  for 

»  .6,  .8,  since  the  entry  is  zero  for  P^  ■  .4).  In  Table  1, 
the  monotonicity  in  r  for  fixed  p  is  the  easiest  to  use. 

Since  f  enters  the  formula  for  temperature  in  a  negative  expo¬ 
nent  (see  formulae  (33)  and  (34)),  the  table  entries  should  improve 
in  accuracy  (for  fixed  p  )  as  T  gets  larger. 


An  Example  of  Table  Use 

To  take  advantage  of  Eq.  (l8),  let  us  note  that 

(p  T,  aQ,  0^,0^, 0^,0^,  c tj)  corresponds  the 

/  (0* )  °C  ..  1  1  \ 

'P«j_l  >  ^  2  *  ) 


to  the  set 
symmetric  set 


Here,  the  index  1  runs  over  0,  1,  2,  3,  4,  3  and  0'  m  2  if 
0-1  and  0'  -  1  if  0  ■  2.  The  correspondence  between  p^'  and 
x-values  is  the  natural  one  (Eq.  (35))>  1 


Suppose,  then,  we  have  the  initial  temperatures  *  1200, 

«  600,  -  200,  UQ  -  100.  Suppose  also  that  we  want  the  tempera¬ 

ture  at  V  m  .1  at  the  right  edge  of  the  first  layer  and  that  aQ  -  0, 

■  10,  ac  ■  1,  ■  .1,  aij  -  l,  ■  .2. 

Here,  pj^  ■  p^  -  1  0  (in  first  string  of  entries),  so  that 

'  -  p^ '  »  0  (in  second  string  of  entries). 
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In  addition, 

t%  «  l,  o£  «  10,  o|  »  0,  a*  ■  1,  aj  ■  .8,  a^'  ■  10,  *  1. 

Equation  (17),  vhich,  in  general,  reads 

“0)V1  +  (D.  •  °o)V2  *  (Ub  -  °o>V 

here,  in  this  example,  is 

(*)  u  -  ioo  -  (noo)v1  +  (5oo)v2  +  (ioo)v5. 

There  remains  only  to  look  up  V1  in  Case  I  tables,  Vg  in  Case  II 

tables,  and  in  Case  II  tables,  under  the  primed  values  of  the 

variables  and  parameters.  These  are  atetttuted  in  (*),  and  we  obtain 
the  value  of  U,  thus: 

In  Table  1,  e.  658,  under  *  .20,  we  read: 

v1(i(1),.i)  -  .8326. 

In  Table  2,  same  slot,  we  read: 

v2(i(1),  .1)  -  .0699. 

In  Table  2,  c.  2052,  under  ■  .8,  we  read: 

V?(0(2),  1.)  -  .0974. 

Substituting  in  (*),  we  get  the  desired  temperature 

U  -  1060.  55. 
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*  initial  temperatures  of  first  and  second  layer 
materials 

U. ,  U_  a  temperature >  as  functions  of  x  and  -7%  in  first  and 
second  layers 

U(p'^, t)  "  temperature  (in  EITHHl  layer)  as  a  function  of  p^^aad  1~ 
U, ,  UQ  »  stagnation  temperatures 

V1  ■  temperature  for  Case  I  (tabulated) 

Vg  ■  temperature  for  Case  II  (tabulated) 

Vj  -  temperature  for  Case  III  (can  be  obtained  from  Vg  table) 

-  temperature  for  Case  IV 
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